Transparent, Non-local, Species-selective Transport in an Optical Superlattice 
Containing Two Interacting Atom Species 
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In an optical superlattice of triple wells, containing two mutually interacting atom species in 
adjacent wells, we show that one species can be transported through the positions of the other 
species, yet avoiding significant overlap and direct interaction. The transfer protocol is optimized 
to be robust against missing atoms of either species in any lattice site, as well as against lattice 
fluctuations. The degree and the duration of the inter-species overlap during passage can be tuned, 
making possible controlled large-scale interaction-induced change of internal states. 
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I. INTRODUCTION 

Optical superlattices offer the means to scale up quan- 
tum operations on individual units of ultracold atoms 
and molecules for parallel processing a multitude simul- 
taneously. Tremendous progress has been made towards 
controlled manipulations in double well lattices [ll-Q al- 
lowing targeted change of internal atomic states and 
controlled merger and separation of pairs of atoms or 
their internal states. In this paper, we propose a new 
level of control achievable in a triple-well superlattice; 
wherein two mutually interacting species can be trans- 
ported through each other with very little (precisely con- 
trolled) overlap at any time. The method dubbed CTAP 
(coherent transport adiabatic passage) operates by prin- 
ciples similar to that of STIRAP (stimulated Raman adi- 
abatic passage, [4]) - an atom in a triple- well (three- 
levels) is transferred between the extreme wells (states) 
avoiding the middle well (state). The phenomenon has 
attracted attention in the context of single species in 
stand-alone quantum wells (B4l0l|. but has been demon- 
strated only with photons in waveguides [ll|. Setting 
aside its novelty value as an interesting quantum ef- 
fect, our goal in this paper is to gauge its true appli- 
cations potential by: (i) applying it to a lattice (Lattice- 
CTAP) with a view to scaling- up and parallel processing 
of atomic operations, (ii) considering two intermingled 
interacting species, with controllable overlap during spa- 
tial manipulation, and (iii) ensuring transfer algorithms 
optimized to work with or without a second species, to 
allow for lattice imperfections and for applications to en- 
tanglement. We map out effects of physical factors and 
their parametric dependencies, and show the feasibility 
of experiments. 

Inasmuch as it involves dynamically controlled coher- 
ent transfer of a quantum state to a very specific target 
state, LCTAP can be viewed in the context of optimal 
control of quantum systems, which has seen much re- 
cent interest in the context time-dependent quantum pro- 
cesses such as loading of Bose-Einstein condensates onto 
optical lattices [l2| and of quantum state transport along 



spin chains [13l - ll5| . As in such phenomena, LCTAP cru- 
cially relies on optimizing the evolution of quantum states 
of atoms in a superlattice, in achieving nonlocal transport 
that avoids obstructions in the classical path. 

The paper is organized as follows: In Sec. [TT] we de- 
scribe our physical model and present our primary results 
demonstrating successful LCTAP evolution in a super- 
lattice. In Sec. IIIII we present our analysis for appro- 
priate optimization of the system parameters to achieve 
LCTAP evolution, specifically highlighting the inherent 
compromise necessary among various desired features. 
We demonstrate feasibility of our results in Sec. IIVI by 
relating the parameters in our simulations to realistic 
physical parameters in current experiments. We conclude 
with comments in Sec. [V] on the assumptions and out- 
look of our model. Details of our numerical simulations 
are provided in Appendix [Al 



II. MODEL AND RESULTS 

We consider an optical superlattice with unit cells of 
triple-wells, that could be implemented with lasers of a 
single wavelength A at three orientations [l6[ with shifted 
frequencies. The lattice potential, 

V(x,t) =4E R [Pi(t)cos[27r(§a; + iyt))] 

+P 2 (t) cos[27r(f a; + P 5 {t))] + P 3 (t) cos(4ttx)] , (1) 

is a superposition of periodicities A/2, 3A/4, 3A/2. The 
lattice depths Pi,P 2 ,P3 (scaled by 4x recoil energy, 
Er = h 2 /2mA 2 of target species A) and the two relative 
phases Pi,P^ (scaled by 27r) constitute five independent 
parameters to be varied in time for LCTAP. We first ex- 
amine and optimize LCTAP for a single target species 
(A) in a superlattice before generalizing for the presence 
of spectator species (B). An approach based on the lowest 
three eigenstates, ground (g), dark (d) and excited (e), of 
each unit triple-well cell, is suitable since adjacent unit 
cells are kept well isolated. LCTAP hinges on maintain- 
ing the system in the dark state @ , since the dark state 
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FIG. 1: (Color online) LCTAP of a single species: Snapshots 
of population density (solid red, left axis) and lattice potential 
(blue line, right axis) for three adjacent super-cells over two 
cycles. The three lowest eigenstates are shown for snapshot 
3 (t=45.8 t). Snapshot times are marked on the evolution of 
the amplitudes Pi,Pi,P3, scaled by 4Er, and relative phases 
P±,P5, scaled by 2n. Movies of evolution are available online 



is antisymmetric with a node in the middle of the cell 
(as can be seen in one of the frames of Fig. [TJ, leading 
to minimal population in the central-well as compared 
to the other two states. Our numerical results assume 
length, energy and time units of \,Er, and t = h/Ej^. 
The time evolution was done by numerically propagating 
the Schrodinger equation for the lattice potential Eq. (TTJ) 
using the split-step operator method based upon stan- 
dard fast Fourier transform (FFT) routines. 



A. Single Species 

In the case of LCTAP with a single species (A), the 
lattice is pattern-loaded with one atom in the left-most 
well of each unit cell [U, [Hj]. Our optimized LCTAP 
achieves simultaneous transfer of all the atoms from well 
1 (left) to well 3 (right) without significant occupation 
of well 2 (middle) during the transfer. We demonstrate 
sustainability by repeating transport from well 3 to well 5 
(well 2 of the next cell) avoiding well 4. Snapshots of the 
transport in Fig. [T] show atom density in three adjacent 
unit cells, and the corresponding lattice potential. Low 
leakage out of each cell and low middle well population 
at all times are evident in plots Fig. |3][a) of evolution of 
the atom density in a cell. 
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FIG. 2: (Color online) LCTAP with two species: Snapshots 
P3] of population density (left axis; broad solid red: species 
A; narrow solid cyan: species B) and lattice potential Va (blue 
line, right axis) shown for one unit cell and one cycle. The 
effect of A-B interaction appears as a 'bump' in the potential 
seen by A; shown separated out for snapshot 3. Evolution of 
lattice parameters Pi is the same with or without the atom B 
present. 



B. Dual Species 

Each unit cell is loaded with an atom A in well 1 and 
an atom B in well 2. Species-selective trapping 0, H3| 
is used to achieve much deeper trap-depth for species B 
(Vb/Va ~ 10), to keep its density profile unaffected by 
lattice variations during LCTAP of species A, which in 
evident in the snapshots in Fig. [2j 

The static shape of the atom B profile justifies a single 
particle approximation, treating the A-B interaction as 
an added external potential felt by atom A. We model it 
on the ID bosonic hard-core interaction giD$(xA — Xb) 
where g\r> — lahuj^ for a tight cylindrical transverse 
harmonic confinement [2l| of frequency oj±, scattering 
length a. Allowing for the density spread, we replace the 
delta function with the ground state density of atom B 
in the middle well, 8{xa ~ xb) — > \ipo(x)\ 2 , so the effect 
of the interaction on A is modeled as a potential bar- 
rier (for positive a) with position and width of \tpQ (x)\ 2 
and strength determined by giD\ipo ( x b)\ 2 - The barrier 
appears as a 'bump' in the potential (Fig. [2]) seen by A. 

The first four snapshots in Fig. [2] demonstrate LCTAP 
of atom A from well 1 to well 3 avoiding overlap with 
atom B in well 2 during the evolution. In the last three 
snapshots atom B is transported oppositely to restore 
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initial configuration (shifted), for the next cycle - also 
achievable by manipulating A alone leaving B in place, 
but requires more complex parameter evolution. The cru- 
cial element of this transfer is illustrated in Figs.Sta) and 
(b), showing similar LCTAP of species A with or without 
the species B actually present, with the same low toler- 
ances on middle well population of A and leakage, and 
most importantly, achieved for both cases with exactly 
the same parameter variation seen in Fig. [2j 

The results discussed here are based on the approxi- 
mation of tight binding of atoms B in the lattice. The 
approximation is reasonable provided that the interac- 
tion energy between atoms A and B is much smaller than 
the energy necessary to excite atom B to a higher vibra- 
tional state. If this approximation is relaxed, the dy- 
namics would contain entanglement of the atomic spatial 
variables during the transport, and the dimensionality of 
the system Hilbert space would increase. This more com- 
plex scenario can have interesting implications that will 
be explored elsewhere, but is not necessary for the ba- 
sic two species LCTAP mechanism we study here, using 
species-specific lattice potential to provide the necessary 
tight confinement of atoms B. 



population at all times (ii) minimal tunneling out of each 
triple well cell (iii) maximal transfer into third well per 
cycle (iv) robustness against noise and (v) maximal trans- 
fer speed. 

Effective LCTAP hinges on the non-adiabaticity factor 
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defined by instantaneous states and their energies. Two 
distinct dependencies are manifest: (i) the energy separa- 
tions AEi = Ei — Ed in the denominator, independent of, 
and (ii) the numerator, linearly dependent on the rate of 
the change of the Hamiltonian. The time rate of change 
can therefore be factored out Aj — £j^f, with s being 
a configuration index labelling the instantaneous states. 
Thus, the optimization is done in two steps: (i) one finds 
the optimal sequence of static instantaneous states and 
(ii) optimizes the dynamic transfer rate factor along that 
path. 



A. Setting Dynamic Conditions 



III. ANALYSIS AND IMPLEMENTATION 

LCTAP is made challenging by the fact that there are 
more physical features to be optimized than the five de- 
grees of freedom offered by the lattice parameters. Here 
we identify the influence of those features and their para- 
metric dependencies, to achieve: (i) minimal middle well 




FIG. 3: (Color online) LCTAP for a single species correspond- 
ing to Fig.[T]shown for one cell and one cycle: (a) Evolution of 
the population density in the three wells; the flat dotted line 
is the leakage out of the unit cell, (b) Evolution of the ener- 
gies with dark state as reference; (c) Evolution of the barrier 
heights and well depths indicated in (d). (e) Evolution of the 
rate-independent adiabatic factors for the ground and excited 
states as a functions of index (s) labeling the instantaneous 
states. 



An intuitive connection to the triple-well case can be 
made by reparameterizing the system in terms of local 
parameters, such as the set of barrier heights and well 
depths {B i: i e [1, 2, 3, 5, 6]} of individual triple- well unit 
cells, shown in Fig. GUd), where the middle well depth, 
£?4, serves as reference. Within our localized state ap- 
proximation, the path for optimal LCTAP evolution can 
be defined precisely in terms of these local parameters. 
Then the corresponding optimal path of the lattice pa- 
rameters {-Pi, i £ [1,5]} is determined by an iterative 
numerical procedure, which is described in detail in Ap- 
pendix [A] The evolution of the local parameters is di- 
rectly dictated by our optimization conditions discussed 
below, and are shown in Fig. (3{c) for the case of sin- 
gle species. The corresponding evolution of the original 
parameters, shown in Fig. [TJ are numerically computed 
from these using the prescription in Appendix [AJ 

By means of this reparameterization, determining the 
optimal LCTAP path translates to finding the optimal 
evolution of an appropriate set of five local parameters, 
that can be mapped to the original five lattice parame- 
ters. As with STIRAP, minimally two time-varying pa- 
rameters are necessary for the transfer. This is accom- 
plished by counter-intuitive variation of the internal bar- 
rier heights f?3 and B§ that fix the couplings of the wells, 
as can be seen in Fig.[3tc), where B$ is lowered before B$. 
The outer barrier height B\ is kept constant in analogy 
with triple wells. 

In order to limit the middle-well population, we im- 
pose two dynamic conditions at every instant: (i) First, 
we keep dark state energy equidistant from the energies 
of the neighboring states A 2 E — AE e — AE g = 0. It 
can be seen from Eq. © that this reduces transitions 
out of the dark state, (ii) Secondly, we maintain the 
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node of the dark state at the center (xq) of the middle 
well, {xo\d} = 0. Geometrical considerations show that 
this lowers the middle-well population due to the anti- 
symmetric shape of the wavefunction at the node. Both 
of these conditions can be implemented with appropriate 
variation of the outer well depths B2 and Be |9| , therefore 
in our simulations, it was convenient to directly reparam- 
eterize the lattice parameters {Pi,i S [1,5]} in terms of 
the parameter set {B 1: B 3 , B 5 , A 2 E, (xo\d)} that define 
the dark state, rather than using the intermediate set 
{Bi,i € [1,2,3,5,6]}. Implementing the two dynamic 
conditions then simply requires setting A 2 E — and 
(xo\d) — for all points on the path. 

B. Compromising on Initial Conditions 

While the two dynamic conditions set the relations be- 
tween the parameters, the actual values are set by initial 
conditions defined by our tolerances for the optimized 
physical features. The key point is to recognize that 
all the undesirable features are maximum at the halfway 
point of a cycle (Fig. [3Jd)) when each cell is symmetric 
about a vertical line at its center. So, parameters chosen 
to satisfy tolerances there will do so over the entire path. 
In this symmetric configuration, three of the five param- 
eters are used to maintain the dynamic conditions: (i) 
dark state energy equidistant from the neighboring en- 
ergies, (ii) dark state node centered in the middle well 
and (iii) the dark state exactly antisymmetric with equal 
population in the extreme wells; the last two correspond 
to setting relative phases P4 = P5 = 0. 

There are two free parameters remaining; however 
there are at least four additional features that we would 
like to optimize: (iv) maximum tolerance on central well 
population, a m (v) the maximum tolerance on the rate 
of leakage out of each cell, (vi) robustness against noise 
and (vii) magnitude of separation AE(= AE e = AE g , 
both being equal in our scheme) between the energy lev- 
els. These four factors are inter-related and obviously 
all cannot be optimized with two parameters; so a right 
balance needs to be found as we now discuss. 



C. Leakage 

One key distinguishing factor that differentiates lattice 
case from CTAP in a stand-alone triple well, is that the 
outer barriers of each unit cell cannot be made arbitrarily 
high due to the interplay of all the lattice parameters in 
defining the barrier heights and well-depths. This leads 
to finite rate of dissipation or leakage out of atom den- 
sity out of each cell into adjacent cells. We find that a 
good measure of the tunneling or leakage out of a cell is 
given by the width of the Bloch band for the dark state 
r d = E d (q = 0) - E d (q = 2tt/3A), where q is the Bloch 
vector for the lattice. This is primarily controlled by set- 
ting the height of the outer barrier B\ between unit cells. 



The relation of Y d with B\ , although monotonically de- 
creasing, is not linear. 



D. Noise Sensitivity 

We cannot predict the exact character of the imperfec- 
tions of superlattices in actual experiments and the pre- 
cise effect they would have on LCTAP evolution, until 
actual experiments are done. But here, we present some 
basic estimates of the tolerance of the process on noise 
and imperfections inevitable in experiments. Fast fluctu- 
ations in lattice parameters Pi can induce non-coherent 
non-adiabatic coupling between eigenstates, which leads 
to a loss of population from the dark state to the neigh- 
boring states. The mean non-adiabaticity (there is no 
sum over i) 

Pl dt 2 ^ (Ej - E d ) 2 dui dt { ' 

j=e,g 

induced by any unwanted fluctuations i/j in each of the 
time- varying lattice parameters Pi is captured by the co- 
efficient /3i and depends on the rate of those fluctuations 
We found that by far, the dominant source of noise 
induced coupling between eigenstates arises from /J4 and 
/?5, associated with the relative phase parameters of the 
lattice potential. Plot of /?4 in Fig. EJa) shows that ef- 
fect of noise is reduced by parameter values that lead to 
higher leakage (T d ) and higher middle well population 
(a m ) both of which are undesirable, thereby indicating a 
compromise in the choice of parameters. 




W Time (units of x) Instant. Config index, s 

FIG. 4: (Color online)LCTAP with two species corresponding 
to Fig. [21 shown for one unit cell and one cycle: Evolution of 
the population density with parameters optimized for pres- 
ence of B in the middle well with atom B (a) present and 
(b) absent, (c) Evolution of the energies of the three states 
with and without atom B (the dark state energies, different 
in magnitude in the two cases is used as reference for both for 
easy comparison) (d) Evolution of the rate-independent static 
adiabatic factors for the ground and excited states, with and 
without atom B, as a function of index labeling the instanta- 
neous states. 
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the product Aj — £j^§ that needs to be within tol- 
erances, LCTAP can be speeded up significantly, by 
keeping Aj fixed and dynamically adjusting the rate 
ds/dt — max(.4 e , „4 g )/max(iS e , S g ); this can double the 
cycle speed, significant in a coherent mechanism. 



FIG. 5: (Color online) (a) Coefficient of noise sensitivity /?4, 
associated with relative phase P4, plotted versus middle well 
population a m and the width of the dark-state Bloch band Yd 
(measure of leakage rate out of unit cell); (b) Plot of inverse 
square of energy gaps between the dark and adjacent states. 
Parameters used in our simulations are marked in both; their 
locations in the middle of the surfaces indicate trade-offs nec- 
essary in optimization. Here Td and AE are in units of Er. 



Our analysis is not dependent on the specific noise 
model. We have chosen 1// noise to probe the influ- 
ence of generic noise in our model, since it is among the 
most ubiquitous in physical systems. An explanation of 
the widespread occurrence of such noise based on the 
self-organized criticality can be found in Ref. (22j. To 
all the parameters P{, we therefore added the frequency- 



dependent noise v oc 



Er=l(/n)" lcOS (27T/^ + . 



with randomly generated frequencies f n and phases c 
which simulates technical noise with spectral density 1//. 
We found that about 3% noise = 0.03 for all i) still 
leads to successful LCTAP, although the transfer profile 
is somewhat distorted. 



E. Energy-level Separation 

Among all the parameters involved in LCTAP, the 
most important one is the energy separation AE between 
the dark state and its adjacent states, since it is related 
to almost all the optimization factors, and Eq. ([2]) sug- 
gests larger AE would be better. But this corresponds 
to lowering the internal barriers B3 and B§ which leads 
to a higher slope of the dark state in the neighborhood 
of the node, leading to higher occupation of the middle 
well. This is a static effect defined by the shape of the 
instantaneous eigenstates, not to be confused with effects 
of transitions induced by nonadiabticity - which do ben- 
efit from a higher energy separation of the states. So a 
compromise has to be struck; Figs. [5] show that if noise 
and energy separation are to be kept optimal as well, we 
need to sacrifice on tunneling and maximum middle well 
population. We chose our parameters by minimizing the 
cost of trade-offs among the optimized features. 



F. Time Evolution 

The simplest possibility for the dynamical factor ds/dt 
is to keep it constant, then the tolerances are de- 
termined by the maximum of the instantaneous fac- 
tors Sj (Fig. GJ|d)). However, recognizing that it is 



G. Dual Species 

All of our analysis applies when species B is present 
in the middle well, but with optimization done for the 
modified potential felt by A. However, to use the same 
LCTAP parameter evolution for both with and with- 
out the atom B present, we change the energy bal- 
ance condition to A 2 E = AE e (without species B) — 
AE g {with species B) = 0, which takes the two most 
proximate levels (as seen in Fig.UJc)) and maintains them 
equally apart from the dark state. The dynamic transfer 
rate ds / dt is now set by the maximum of four adiabatic 
factors (Fig. HJd)) corresponding to the ground and ex- 
cited states, with and without the atoms B. 



IV. PHYSICAL PARAMETERS AND 
FEASIBILITY 

We now put our parameters in the context of recent 
experiments with 87 Rb in superlattices, choosing laser 
wavelengths around 803 nm close to the Dl and D2 lines. 
This translates to time unit r = 45 microseconds putting 
our LCTAP cycles at 2.5—4.5 milliseconds near accessible 
time scales [2j ■ Figures Q] and [2] show that the individual 
lattice depths (Pi,P 2 ,P3 in units of 4_Er) for the target 
species A range between 2E R — 8E R , their sum never ex- 
ceeding about 10E R . We found that to keep the species 
B unaffected as target species A evolves, Vb/Va — 6 
is sufficient [l7j implying a maximum net lattice depth 
for species B of about 60E R , accessible in current ex- 
periments 0, 0|. Figures [3] and 0] use a higher value 
Vb/Va = 20 only to show the interaction potential as 
an obvious 'bump'. 

The ground state density of atom B is a Gaussian of 
a = 0.035A and peak density IV^O)! 2 = n -4/A- For our 
estimates, we assume two hyperfine states [2fJ| of 87 Rb 
and use its triplet scattering length a = 99a b- For a 
uniform transverse lattice confinement equal to Vb — 
60E R we get w_l = (2n/X) v / 2V B /m = 3.4 x 10 5 Hz, 
gin = 0.20E R X and barrier height g 1D \ip^ {0)\ 2 = 2.3E R . 
It can be much less for a transverse magnetic confinement 
on an atom chip [23[, wj_ ~ 2?t x 5100 Hz and gip = 
0.019E R \. Our simulations use a cautious higher value 
for the barrier height 0.30^(0) | 2 = 3.4£^. 

Finally, we note that high transfer efficiency was 
achievable in our simulations. In the case of single 
species, shown in Figs. Q] and [3l at the end of the first 
cycle 94.6% of the population ended up in the right well 
with 4.7% overall lose out of each cell, the remainder 
distributed between the first and the middle well. With 
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dual species, as shown in Figs. [2] and |U the correspond- 
ing numbers were 90.8% and 6.3% with species B present 
in the central well; and 91.3% and 7.5% with species B 
absent in the central well. These values indicate that 
multiple cycles can be repeated in a single run of an ex- 
periment. 

V. CONCLUSION 

We note that our use of a three-mode localized state 
is validated by our numerical results: (i) there is low 
leakage of density out of each cell per cycle and (ii) the 
width of the dark state Bloch band T d ~ 0.01E R (Fig.HJ 
is an order of magnitude less than the energy separations 
(Fig. [3] and 2]). Our results indicate viable implementa- 
tion of LCTAP and should motivate experiments with 
triple-well optical lattices. By allowing choice of the de- 
gree and the time of overlap of the two species this mech- 
anism can be tailored for controlled interaction-induced 
change of internal states. 
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Appendix A: Numerical reparametrization 

A key step in our simulations is the reparameterization 
of the lattice variables {Pi,i € [1,5]} of amplitude and 
phase in terms of characteristics of the localized states of 
individual triple-well unit cells. 



The lattice Hamiltonian H = H(Pi(t)) is characterized 
by five time-dependent lattice parameters, {Pi, i £ [1,5]}. 
We identify five characteristics {Fj,j £ [1,5]} of the 
Hamiltonian that need to be optimized for LCTAP. The 
Hamiltonian can be numerically reparametrized in terms 
of these properties Fj(t) = F J [H{P t {t))] = Fj[P z (t)] pro- 
vided that the Jacobian matrix Jji = has nonzero 
determinant in the region of the reparametrization. 

The properties Fj are our operational parameters, 
with the advantage that our optimization conditions for 
LCTAP directly determine the dynamic values of these 
new parameters Fj. In our numerical simulations, where 
the evolution is done by split step operator method, at 
every time step t n of the evolution we therefore have a 
well defined set of known optimal parameters Fj(t n ). In 
order to determine how the lattice needs to evolve for op- 
timal LCTAP, at each of these steps we need to determine 
the corresponding optimal lattice parameters Pi (t n ) . We 
do this iteratively as follows: 

At every time step t n , start with a reasonable guess 
for the lattice parameters based upon the values of 
the previous time step, but which correspond to non- 
optimal values Fj — Fj(P^ ). Using the differences 
SFj — Fj(t n ) — Fj in the system of linear equations 

«f } = E l = Ei JjM 1] , we solve for 5P™ 

by inverting the Jacobian matrix, and thereby find a new 

improved set of parameters P^ = P^ + 5P^ with the 

(2) (2) 

new values Fj ; = Fj (P> ; ) closer to the optimal values 
Fj(t n ). We repeat this for k iterations until SF^ = 
Fj(t n ) — Fj is less than our set tolerance. 

As described in Sec. IIII1 the appropriate properties Fj 
are chosen to be {B\, B3, B$, A 2 E, (xo\d)}. As discussed 
in Sec. IIII1 we set A 2 E = and (xo\d) = for all points 
of the path to meet our conditions for LCTAP. 
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